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In this paper, we investigate a family of one-dimensional multi-component quan- 
tum many-body systems. The interaction is an exchange interaction based on the 
familiar family of integrable systems which includes the inverse square potential. 
We show these systems to be integrable, and exploit this integrability to completely 
determine the spectrum including degeneracy, and thus the thermodynamics. The 
periodic inverse square case is worked out explicitly. Next, we show that in the limit 
of strong interaction the "spin" degrees of freedom decouple. Taking this limit for 
our example, we obtain a complete solution to a lattice system introduced recently 
by Shastry, and Haldane; our solution reproduces the numerical results. Finally, we 
emphasize the simple explanation for the high multiplicities found in this model. 



The plan of this paper is to first review integrable systems and the Lax technique, show 
that certain modified systems with an exchange interaction are integrable, give the general 
solution to these systems and an explicit solution for an important example, then to produce 
from this solution by an appropriate limit the solution to a family of lattice problems, and 
finally we discuss what it all means. The impatient reader may skip to the closing paragraphs 
for a summary of the results. 

I. INTEGRALS OF MOTION 

If we consider a classical system of N one-dimensional particles governed by the Hamil- 
tonian 

H = \ P)+ J2 vixj-Xk), 

l<j<N l<j<k<N 

then it has been shown by Lax jl], Moser || and Calogero |§ that for certain potentials one 
can find two Hermitean N x N matrices L and A that obey the Lax equation 

dL/dt = i[A,L\. 

Thus L evolves by a unitary transformation generated by A, and hence det[L — wl] is a 
constant of motion. Expanding the determinant in powers of uj we find N integrals of 
motion 

det[L-uI}= Jjh^) N - j , 

0<j<N 

j = l,...,N. Further, these integrals have been shown to be in involution, and thus the system 
is integrable. 

The two matrices are given as 

Ajk = 5 jk ^ l( x j ~ x i) + (1 - $jk)(3(xj ~ Xk), 
Ljk = SjkPj + i(l - 5 jk )a(xj - x k ). 

The functions v(x), 7(x) and /3(x) are even, while a(x) is odd. They obey the equations, 
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v(x) = —2a(x)fl(x), 
/3(x) = —a(x), 
a(x + y)[y(y) - 7(2;)] = a(y)a(x) - a(x)a(y). 

Calogero || has shown that the most general solution to these equations is given in terms 
of the Jacobi elliptic function sn(x\m) as 

a(x\m, A, k) = Xa(nx) 



sn(K,x\m) ' 



Xk 

j(x\m, A, k, c) = Xk^(k,x) + c = Xna 2 (nx) + c = — — : — ; — - + c, 

sn 2 {Kx\m) 

v(x) = X 2 a 2 (K,x) + v . 

Calogero || has also demonstrated that if one replaces the classical dynamical variables with 
the corresponding quantum mechanical operators, det[L— wl] is well defined with no ordering 
ambiguity, and the quantum mechanical commutator \H, det[-L — wl\\ = 0. Thus, the Jj 
are still constants of motion. Finally, Calogero showed that [ det[L — wl], det[L — w'l] | = 0, 
and thus the quantum system is also completely integrable. 

To sum up the situation: The one dimensional quantum system, governed by the Hamil- 
tonian 

is completely integrable. For the general classical system integrability tells us something 
concrete, namely that the motion in terms of action-angle variables is on a torus. However, 
for the general quantum system integrability seems to buy one almost nothing. The excep- 
tion is for those special cases which support scattering; i.e., systems which fly apart when 
the walls of the box are removed. In these cases, in the distant past and future the Lax 
matrix L approaches a diagonal matrix, so that 

det[L — col] = Yl (Pj ~ UJ )- 

Thus the individual momenta pj are conserved in a collision, and hence the wavefunction is 
given asymptotically by Bethe's ansatz. Sutherland || has exploited this fact to completely 



determine, in the thermodynamic limit, properties of systems interacting by potentials of 
the forms v(x) = g/x 2 ,g/ sin 2 (a;) and gj sinh 2 (a;), including the Toda lattice. We emphasize 
that no features of the proof of integrability are needed. All we need is to know it to be 
integrable, by whatever method. (We might even just assume it to be integrable, and deduce 
the consequences.) 



II. A MODIFIED PROBLEM 

Although the previous problems are multi-component problems, in fact statistics enters 
only trivially, due to the strong repulsion of the potential at the origin. Polychronakos 



has modified this problem |]10Hll| , allowing particles to penetrate by means of an exchange 
interaction. He then proves the integrability of this modified problem. We now offer an 
alternative proof of integrability; we extend and elaborate on this method in a separate 
paper [|12|]. If Pj k is the permutation operator that permutes particles j and k, let us modify 
the Lax matrices to be 

A jk = 5 jk Pjil( x j ~ x i) + (1 - Sjk)Pjkf3(xj - x k ), (1) 
Ljk = SjkPj + i(l - 5 jk )P jk a(xj - x k ). (2) 

We assume that the new potential Vj k commutes with Pj k . We then seek to satisfy a 
quantum Lax equation of the form [H, L[ = [L, A]. Here the first fancy commutator is a 
quantum mechanical commutator between operators, while the second commutator is an 
ordered matrix commutator, so that the equation above is really N 2 equations of the form 

[H,Lj k }= {LjiAi k — AjiLik). 

1<1<N 

In order that these equations be satisfied, one finds that the functions a{x),(3{x) and 7(2) 
must satisfy exactly the same equations as before. The potential now is 



Vj k (x) = X a (kx) + \K^(Kx)Pj k 



sn 2 (/tx|m) 
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We will henceforth be interested in a particular case. Let k — i, X — i, c — — 1 and 
m — > 1, so that 

a(x) — > cot(x), 

7(2;) — > —1/ sin 2 (x) = — 1 — a 2 (:r), 
a' (a;) = —@(x) — > —1/ sin 2 (x) = 7(2;). 

Restoring scale factors, we have a potential 

Vjk( x ) = (A 2 — XPjk)^ 2 1 sinh 2 (/ta;). 

(We could as well have a sine instead of a sinh.) 

The significance of this case is for the form of the Lax A matrix. Since (3{x) = —7(2;), 
then defining a vector r\ with 77^ = 1, we see Aq = rfA = 0. This allows us to construct 
constants of motion by I n = L n r), since 

|lf,J B |=i7t|lf,L»to = i7t ]T L|L JV ~ 1- ' J }?7 

0<i<Af-l 

= 77+ 53 {V[A, L]L N - x -i}ri = if {AL N ~ l - L N ~ 1 A}r] = 0. 

0<j<AT-l 

By Jacobi's relation for commutators, |7 n ,7 m J is a constant of motion, and since this is a 
system that supports scattering, we see \I n , I m \ — > 0, and hence the system is completely in- 
tegrable. Thus once again, the individual momenta pj are conserved in a collision, and hence 
the wavefunction is given asymptotically by Bethe's ansatz. This allows one to completely 
determine, in the thermodynamic limit, properties of systems interacting by potentials 

v (x) = (A 2 - XP)/x 2 

and 

(A 2 - AP)/sinh 2 (:r), 
as well as the periodic versions of these such as 

(A 2 - XP)-n 2 /L 2 sin 2 (7nr/L) 
in the limit as L — > 00, density finite. 
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III. THE (A 2 - XP)/X 2 AND (A 2 - XP)tt 2 /L 2 sin 2 (7rX/L) PROBLEMS 



The input for the Bethe ansatz wavefunction is the two-body scattering matrix, since the 
A-body scattering is simply a product of all N(N — l)/2 scatterings. At the same time, we 
can gain some confidence by confirming the consistency of the two-body scattering matrix; 
it must satisy the Yang-Baxter equations. With the relative coordinate r = x 2 — Xi, the 
potential is 

v (r) = (A 2 - AP)/r 2 , 

so depending on whether the wavefunction is even or odd, the potential is v±(r) = A(A =F 
l)/r 2 . The radial equation then is —</?" + a(a — l)ip/r 2 = k 2 ip, with a = A or A + l according 
to whether the wavefunction is even or odd. The solution is given as 

(p(r) = y/rJ a - 1/2 (kr), r > 0, 
— > cos(A;r — ira/2), r — > oo 
— > r a , r — > 0. 

For the even and odd solutions then we have y?±(r) — > e~* fcr ± e -* 7rA e ifer ? r _> oo. 

To construct a scattering state in which particle 1 with momentum k\ is incident on 
particle 2 with momentum fc 2 , ki > k 2 , we take the linear combination 

_^ e -i W Agi(fclXl+fc2X2) j X2 « Xl 

Thus, there is no reflection, and the particles pass through each other with only a phase 
shift e _l7rA . Consistency is not an issue. This is a surprising result with important conse- 
quences. 

We now place the particles on a large ring of circumference L, making the problem 
periodic. We thus need to replace the potential with some periodic potential, such as 
v(x) = (A 2 — \P)ir 2 /L 2 sm 2 (nx/L), which reduces to (A 2 — \P)/x 2 in the limit L — > oo, and 
which therefore has the same thermodynamic limit. To repeat: We will be giving the exact 
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solution for the inverse sine squared exchange potential in the thermodynamic limit of an 
infinite system with finite density. 

Taking the jth particle around the ring of circumference L, it scatters from every other 
particle with a net phase shift that must be unity, to insure that the wave function is periodic. 

Thus the k's satisfy the equation exp[ikjL — inX s g n (^j — h)] = 1, or upon taking the 

'(/./) 

logarithm, kj = k° + ^ J2 sgn(kj — k{). 

Here k° = 2tt ■ integer • L are the the non-interacting momenta with k° > k% > ... > k° N . 
We easily solve for the fc's obtaining kj = k° + 7rA(iV + 1 — 2j)/L. Then the energy is 

E = \ £ k) = \ £ (k°) 2 + ^ £ k°(N+l-2j) + ^N(N*-l)/QL^ 

1<3<N l<j<N ^ l<j<N 

while the momentum is 

P= E *; = E k °r 

l<j<N l<j<N 

All of this is just a repetition of the old results H for the usual 1/r 2 potential, which in 
fact is a special case of this interaction with a single species of particle. 

Although the energy eigenvalues are exactly same, these levels are now highly degenerate. 
In fact, the degeneracy are exactly the same as the free particle degeneracy when A — > 0. 
Thus, let us introduce occupation numbers v a (k), equal to the number of particles of species 
a with non- interacting momenta k. Thus v a {k) = 0, 1 for fermions, and v a {k) = 0, 1, 2, . . . 
for bosons. Then the energy E is given by 

j /*oo \ /*oo fOO 

E/L = —y2 dk v a {k)k 2 + — V / dkv a {k) dk' v p {k')\k- k'\ + 7c 2 \ 2 d 3 /6 

47T ^ J- oo o7T J —oo J — ■ oo 

1 roo \ roo roo 

= — dkv{k)k 2 + — dku(k) dk' u(k')\k-k'\+7T 2 \ 2 d 3 /6 

4:71 J-oo 07T J-oo J-oo 

Here, the density of species a is d a = N a /L = ^- J^dk v a {k) , v{k) = ^v a {k) , and 

a 

d = N/L = d a — y~ I^oo d& v(k) is the total density. The degeneracies are now expressed 

a 

in terms of the usual entropy S, as 

roo 

S/L = (1/2tt) £ / dfc (^X 1 ± "«) ln (! ± v a) ~ "a ln(z/ Q ). 
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(Here, and in what follows, the upper sign is for bosons, while the lower sign is for fermions.) 

We can now determine the equilibrium occupation numbers v a {k) as those which maxi- 
mize 

P/T = S/L - E/TL + J2 VaN a /TL. 

a 

This gives the familiar expression is a = ]_/^ e ~^ a ^ T =F 1]. In this expression, 

/oo 
dib' v(k')\k-k'\+7r 2 \ 2 d 2 /2. 
-oo 

Differentiating, we find for the interacting momenta p that appear in the asymp- 
totic Bethe ansatz, p = e' = k + (A/2) J^dk' u(k')sga[k — k'\ . Differentiating once 
again, p' = e" = 1 + \v = p& = ^ . Thus, integrating, we obtain p 2 /2 = 
E + XTJ2 (±1) ln[l =f e^ a_£ ^ T ]. We may then go back and write the densities as 2nd a = 

a 

f™oc dpjf 1 ^- This gives us a complete solution for the thermodynamics. 

For instance, if we consider a system consisting of two species of bosons, say + and 
-, then taking /x± = fi ± h, we find e ~ p2/2XT = [e^-^^ + cosh(/i/T)] 2 - smh 2 (h/T), so 
e (e-n)/T = { e ~ P 2 /2\T + sinh^/i/T)] 1 / 2 - cosh(/i/T) . This then allows us to write u± as an 
explicit function of p, which may be integrated to give d± as a function of h±, and which in 
turn can be integrated with respect to h± to give the free energy, which can be referenced 
to the known one-component system at h + — > oo. 

IV. GROUND STATE ENERGY 

At this point, let us look at the ground state energy as a function of the densities 
{d a }. This is easily done, and it is only for the ground state when T = that we expect 
singularities. We assume / species of fermions with densities d\ < c?2 < • • • < df, and a 
total of d B bosons of any species, so the total density is d. Then the occupation numbers 
are the free particle ground state occupation numbers, and the integrals in the expression 
for the energy are easily evaluated to give 

QE/ir 2 L = (1 + 2A) ]T d 3 a + 3\d B £ d 2 a + X ]T d a (d 2 a + 3d}) + \ 2 d 3 . 

l<a<f l<a<f l<Q</3</ 
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Because this expression is only valid for a particular ordering of the densities, the ground 
state does have singularities whenever d a = dp. A Taylor expansion gives the singularity 
as E/L « \d/3 — d a \ 3 , so the susceptibilities do not diverge. The singularity structure, and 
even the nature of the singularities, is very similar that found previously for a short-ranged 
lattice exchange model J7J. 

V. THE LATTICE PROBLEM 

Examining the potential energy as A — > oo, we see that the particles crystallize into a lat- 
tice with lattice constant 1/d, and so the elastic modes of this lattice, and the compositional 
or "spin" modes on the lattice separate. We thus have as A — > oo, 

H (1/2) £ pj + A* £ — ., fl -Aj £ >ra ^.,r' P,> 



H e + Xd 2 H s . 



i<jTk<N L2 sin^Tr^- - x fe )/L] i<^<jv ^ 2 sin 2 [7r(j - fc)/JV] 



The Hamiltonian for the elastic modes is just the familiar one-component system ||, al- 
though with a coupling constant A 2 . We have defined the spin Hamiltonian H s with ex- 
change constant unity, although the effective exchange constant of the last term of H is \d 2 . 
The two-species problem was first studied independently and simultaneously by Shastry ||, 
and by Haldane M. Thus we see that by an appropriate limiting process we may determine 
the thermodynamics of this spin Hamiltonian exactly. For the spin problem, the appropriate 
concentrations are m a = N a /N = d a /d. 

Let F(T, {N a }, L) be the free energy for the original continuum exchange problem. Then 
it is seen on dimensional grounds to be of the form LT 3 / 2 / ({d a } /T 1 / 2 ) . Given a set of "chem- 
ical potentials" h a , such that h a = , then we make a partial Legendre transformation, 

a 

F(T, {N a }, L)-J2 KN a = F(T, {h a }, N, L) = LT 3 ' 2 f(d/T^ 2 , {K}/T^ 2 ). 

a 

(Free energies are identified by their variables.) Then N a = —dF/dh a , so if d a (d, {h a }) = 
—df(d,{h a })/dh a , we can write d a = T x i 2 d a {d/T x l 2 , {h a }/T), which for T = 1 is simply 
d a . 



Returning now to the limit as A — > oo, and setting T = 1, we have f{{d a }) — > + 
df s (l/Xd 2 , {m a }). Here is the free energy density for the one-component system || 

with coupling constant A 2 , and f s (T, {m a }) is the free energy per site for the spin system. 
Again making the partial Legendre transformation, we find 



Remember, in the functions on the right-hand side, the first arguement represents the 
density d = 1/y/XT — > as A — > oo. Finally, the relationship is a little clearer in terms 
of the concentrations, m a = lim d a /d({fi a }) , where the chemical potentials are {fi a } = 



{/i + h a /T}, and \i is adjusted to give 1 = TXd 2 . 

Referring to our previous explicit results, in the limit as A — > oo, with 1 = TXd 2 , then 
m a = d a /d = 1/27T = f* n dpv a /v, with u a = 1/ [exp(e — h a /T) =f 1] , and {p 2 — tt 2 )/2 = 
TJ2 (±1) m [l T exp(h a /T — e)]. This gives us a complete solution, since the concentrations 



can be integrated with respect to {h a } to give the free energy. The free energy is referenced 
with respect to the trivial case of m a = 1 for some a. For a two species problem of either both 
fermions or both bosons, our results reproduce Haldane's conjecture based on his numerical 



from ferromagnetic to antiferromagnetic simply by changing the statistics of the particles. 



Let us now summarize what we have accomplished. We have investigated a family of 
one-dimensional multi-component quantum many-body systems, interacting by an exchange 
interaction based on the familiar family of integrable systems which includes the inverse 
square potential. These systems are shown to be integrable. Since these systems support 
scattering, the asymptotic wave function is given by Bethe's ansatz, and consequently the 
spectrum can be completely determined, in the thermodynamic limit. The interactions that 
can be solved in this way include v(x) = (A 2 - XP)/x 2 and (A 2 - AP)/sinh 2 (x) at finite 




A— >oo 



a 



investigation [^]. We emphasize that one changes the sign of the of the spin Hamiltonian 



VI. SUMMARY AND CONCLUSION 
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density, as well as the periodic versions of these such as (A 2 — \P)it 2 /L 2 sm 2 (Tix/L) in the 
limit as L — > oo, density finite. 

We explicitly work out the case of the inverse square or inverse sine square potentials, 
and demonstrate that for this exchange potential particles scatter without reflection. They 
pass through one another with only a phase shift, and this phase shift is independent of 
the identities of the particles (and the momentum as well). As a consequence of this per- 
fect transmission, there can be no Bragg reflection to split the energy of the even and odd 
solutions when the particles pass around the ring. Instead we have degenerate states, and 
so we could take a linear combination which, for a two-component system, would have all 
particles of one type going clockwise around the ring, while all particles of the other type 
go counter-clockwise. Yet this state would be an eigenstate of energy. This then explains 
the very high degeneracy observed in the spectrum; it is just the degeneracy of the free 
particle system. We can offer other examples of systems without reflection, with similar 
behavior. In fact, it is clear from the discussion, that by using adiabatic continuity, the 
quantum numbers of the states can be taken as the corresponding quantum numbers of 
the free particle states. Singularities occur only in the ground state, and these are exhib- 
ited. Finally, it is demonstrated that for each of these systems, taking the limit of infinite 
coupling constant appropriately gives the solution to a corresponding lattice problem. One 
could easily prove the integrability of these lattice systems directly using the Lax technique. 
However, application of the asymptotic Bethe ansatz directly to the lattice is plagued by 
umklapp problems; the momenta are pushed by the interaction outside the Brillouin zone 
|10| . The limiting trick avoids these complications. Taking this limit in our explicit ex- 
ample, we obtain a complete solution to a corresponding lattice system introduced recently 
by Haldane, and by Shastry. Our solution reproduces Haldane's numerical results. This 
system exhibits very high multiplicities; we find these multiplicities to be just the familiar 
degeneracy of free particles. It arises naturally from the absence of reflection, and hence of 
Bragg reflection and the consequent splitting of levels, and is nothing more. This seems to 
us a good explanation. 
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The success of the asymptotic Bethe ansatz has been rather dramatic, so let us here 
reiterate the philosophy behind the method. To justify using the asymptotic Bethe ansatz, 
one needs all of the following: a system which supports scattering, a proof of integrability 
with integrals that depend only on the asymptotic momenta, and a valid virial expansion. 
Given these, and using only the two-body scattering data as input, one can deduce all ther- 
modynamic properties, including ground state properties and low-lying excitations, within 
the phase for which the virial expansion is valid. Often other phases may be reached by 
symmetry. 

We conclude by reemphasizing that there are many other interesting systems to be 
explored in detail in the future. 
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